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ABSTRACT 

We investigate the conjecture by Sikora, Stawarz & Lasota (2007) that the observed AGN-radio- 
loudncss bimodality can be explained by the morphology - related bimodality of black-hole spin 
distribution in the centers of galaxies: central black holes in giant elliptical galaxies may have (on 
average) much larger spins than black holes in spiral/disc galaxies. We study how accretion from a 
warped disc influences the evolution of black hole spins and conclude that within the cosmological 
framework, where the most massive BHs have grown in mass via merger driven accretion, one indeed 
expects most supermassive black holes in elliptical galaxies to have on average higher spin than black 
holes in spiral galaxies, where random, small accretion episodes (e.g. tidally disrupted stars, accretion 
of molecular clouds) might have played a more important role. 

Subject headings: cosmology: theory - black holes - galaxies: evolution - quasars: general 



1. INTRODUCTION 

It has been known for many years that the radio loud- 
ness of AGN hosted by disc galaxies is on average three 
orders of magnitude lower than the radio loudness of 
AGN hosted by giant ellipticals (see Xu et al. 1999 and 
references therein). However, as as shown by HST obser- 
vations, such a galaxy morphology - radio- loudness cor- 
respondence is not "one-to-one" : both radio-quiet and 
radio-loud very luminous quasars are hosted by giant el- 
lipticals (Floyd et al. 2004). On the other hand the 
popular version of the so-called "spin paradigm" asserts 
that powerful relativistic jets are produced in AGN with 
fast rotating black holes (Blandford 1990), implying that 
BHs rotate slowly in radio-quiet quasars, which repre- 
sent the majority of quasars. However, such conjecture, 
at least in its basic interpretation, is in conflict with the 
high average BH spin in quasars deduced from the high 
average radiation efficiency of quasars using the "Soltan 
argument" (Soltan 1982; Wang et al. 2006 and references 
therein). 

Parallel studies of radio-emission from X-ray binaries 
showed that at high accretion rates production of jets is 
intermittent (Gallo et al. 2003) and that this intermit- 
tency can be related to transitions between two different 
accretion modes (Livio et al. 2003). This inspired Ul- 
vestad & Ho (2001), Merloni et al. (2003), Nipoti et 
al. (2005), and Hording et al. (2006) to postulate the 
existence of a similar intermittency of jet production in 
quasars and formulate an "accretion paradigm" accord- 
ing to which the radio-loudness is entirely related to the 
states of accretion discs. However, Sikora et al. (2007) 
found that on the radio-loudness - Eddington-ratio plane 
AGN form two parallel sequences that occurrence of can- 
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not be explained by the "accretion-paradigm" (see also 
Terashima & Wilson 2003; Chiaberge et al. 2005; and 
Panessa et al. 2007). Sikora et al. (2007) therefore pro- 
posed a revised version of the "spin paradigm" , suggest- 
ing that giant elliptical galaxies host, on average, black 
holes with spins larger than those hosted by spiral/disc 
galaxies. 

This morphology-related radio dichotomy breaks down 
at high accretion rates where the dominant fraction of 
luminous quasars hosted by elliptical galaxies is radio 
quiet. This radio-quietness occurs in quasars with high 
spin values. In such systems with high accretion rates the 
intermittency is related to the conditions of production 
of collimated jets, in agreement with what is found in X- 
ray binaries, and with the "Soltan argument". It should 
be emphasized that even if the production of powerful 
relativistic jets is conditioned by the presence of fast ro- 
tating BHs, it also depends on the accretion rate and on 
the presence of disc MHD winds required to provide the 
initial collimation of the central Poynting flux dominated 
outflow. 

In this article we will examine under which condition 
the cosmological evolution of BHs in galaxies may lead to 
low spins in disc galaxies and high spins in more massive 
ellipticals. 

To put our investigation in the relevant context we will 
first recall why the value of black-hole's spin might be 
of fundamental importance for relativistic jet launching. 
Assuming that relativistic jets are powered by rotating 
black holes through the Blandford-Znajek mechanism, 
Blandford (1990) suggested that the efficiency of jet pro- 
duction is determined by the dimentionless black hole 
spin, a = Jh/Jmax = cJh/GM^, where J h is the angu- 
lar momentum of the black hole. If true, this could ex- 
i>iap.]&iain the very wide range of radio-loudness of AGN that 
look very similar in many other aspects by attributing it 
to a corresponding black-hole spin distribution. This so 
called "spin paradigm" was explored by Wilson & Col- 
bert (1995), who assumed that the black-hole spin evo- 
lution is determined mainly by mergers. They claimed 
that mergers of black holes, following mergers of galaxies, 
lead to a broad, 'bottom-heavy' distribution of the spin, 
consistent with a distribution of quasar radio-loudness. 



2 



However, this claim was challenged by Hughes & Bland- 
ford (2003), who showed that mergers cannot produce 
the required fraction of black holes with high spins and 
concluded that accretion of matter is essential in deter- 
mining black-hole spins. In this case, however, as no- 
ticed earlier by Moderski & Sikora (1996) and Moderski, 
Sikora & Lasota (1998; hereafter MSL), one encounters 
the difficulty of maintaining a sufficient number of black 
holes at the required low spin, the spin-up by accretion 
discs being so efficient. MSL could match the distribu- 
tion of radio-loudness with the spin distribution only by 
feeding holes with very small randomly oriented accre- 
tion events, i.e. by accretion events forming co-rotating 
and counter-rotating discs with the same probability. 

MSL also addressed the problem of the spin overflip- 
ping due to the Bardeen-Petterson effect. When an ac- 
cretion disc does not lie in the equatorial plane of the 
BH, that is, when the angular momentum of the accre- 
tion disc is misaligned with respect to the direction of 
Jh, the dragging of inertial frame causes a precession 
that twists the disc plane due to the coupling of Jh with 
the angular momentum of matter in the disc. The torque 
tends to align the angular momentum of the matter in 
the disc with that of the black hole, causing thus the in- 
clination angle between the angular momentum vectors 
to decrease with decreasing distance from the BH, forc- 
ing the inner parts of the accretion disc to rotate in the 
equatorial plane of the BH (Bardeen & Petterson 1975). 
Sustained accretion from a twisted disc would align the 
BH spin (and the innermost equatorial disc) with the an- 
gular momentum vector of the disc at large radii (Scheuer 
& Feilcr 1996). If the disc was initially counter-rotating 
with respect to the BH, a complete overflip would even- 
tually occur, and then accretion of co-rotating material 
would act to spin up the BH (Bardeen 1970). 

MSL concluded that the Bardeen-Petterson effect can 
be neglected because the alignment time (10 7 years; Rees 
1978) is longer than the duration of a single accretion 
event. Later, however, a series of papers put into doubt 
the validity of Rees's estimate (Scheuer & Feiler 1996, 
Natarayan & Pringle 1998). This framework was re- 
cently investigated by Volontcri ct al. (2005) who ar- 
gue that the lifetime of quasars is long enough that an- 
gular momentum coupling between black holes and ac- 
cretion discs through the Bardeen-Petterson effect effec- 
tively forces the innermost region of accretion discs to 
align with black- hole spins (possibly through spin flips), 
and hence all AGN black-holes should have large spins. 

Recently King et al. (2005) pointed out that un- 
der some conditions the alignment torque can lead to 
disc-hole counter- alignment reactivating the debate. The 
counter- alignment process was numerically simulated by 
Lodato & Pringle (2006). 

We here re-analyze the alignment problem in view of all 
these latest results. We explore what are the likely out- 
comes of accretion episodes that grow black holes along 
the cosmic history, and determine under which conditions 
black holes in disc galaxies end up having low spins. 

2. ASSUMPTIONS 

The dynamics involving a spinning black hole accret- 
ing matter from a thin disc whose angular momentum is 
not aligned with the spin axis has been studied in a num- 
ber of papers (e.g., Papaloizou & Pringle 1983; Pringle 



1992; Scheuer & Feiler 1996; Natarajan & Pringle 1998). 
A misaligned disc is subject to the Lense-Thirring pre- 
cession, which tends to align the inner parts of the disc 
with the the angular momentum of the black hole. The 
outer regions of the disc are initially inclined with respect 
to the hole's axis, with a transition between alignment 
and misalignment occurring in between at the so-called 
"warp radius" (see below). The direction of the angular 
momentum of the infalling material changes direction as 
it passes through the warp. In our calculations we will 
assume that the black-hole spin evolution is determined 
by accretion only. Volonteri et al. (2005) have shown 
that BH mergers play a sub-dominant role in the global 
spin evolution. 

2.1. Viscosities 

Despite many efforts the problem of warped discs, es- 
pecially in the non-linear regime, has yet to be solved. 
Therefore the characteristic scales of the problem arc 
subject to several uncertainties. The main quantity 
of interest is the "warp radius" R w defined as the ra- 
dius at which the timescale for radial diffusion of the 
warp is comparable to the local dragging-of-inertial frame 
( "Lense-Thirring" in the weak-field approximation) time 
{acR^/R 3 )- 1 (Wilkins 1972). The timescale for the 
warp radial diffusion can be written as 



where V2 is a viscosity characterizing the warp propa- 
gation which can be different from the accretion driv- 
ing viscosity, v\, which is responsible for the transfer of 
the component of the angular momentum parallel to the 
spin of the disc. The relation between v\ and vi is the 
main uncertainty of the problem, assuming of course that 
such two-viscosity description is adequate at all. De- 
scribing v\ by the Shakura-Sunyaev parameter a one 
can show (Papaloizou & Pringle 1983) that the regime 
in which H/R < a <C 1 (H being the disc thickness) 
one has ^1/^2 ~ ct 2 . In such a case the accretion time 
£ aC c ~ Rwjv\ would be much longer than the warp dif- 
fusion time t w . 

However, such a description can be questioned on sev- 
eral grounds. First, is the a < 1 appropriate for high- 
rate accretion onto AGN black-holes? There are no reli- 
able estimates of this parameter for AGN but outbursts 
of LMXBs suggest that in hot accretion discs a > 0.1 
(see e.g. Dubus et al. 2001). In such a case v\ is com- 
parable to i> 2 (Kumar & Pringle 1985). Second, even if 
the two viscosities are different is t w the relevant time 
for black-hole re-alignment? This is not clear since this 
latter process is very dissipative and could be controlled 
by accretion and not warp propagation. 

2.2. Relevant radii 

During the accretion process, the angular momentum 
of the disc at the warp location sums up with that of 
the black hole, so the angle between the angular momen- 
tum of the outer disc and the BH spin changes. King 
et al. (2005) suggest that the condition of alignment or 
counter- alignment can be expressed as a function of the 
angular momenta of the hole and of the disc: Jh and 
Jd- The counter-alignment condition depends on the 
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ratio 0.5 Jd/ Jh, to be compared with the cosine of the 
inclination angle, <j>. If cos</> < — 0.5 Jd/ Jh, the counter- 
alignment condition is satisfied. King et al. (2005), how- 
ever, leave the definition of Jd vague, indeed they sug- 
gest that Jd is the angular momentum of the disc inside 
a certain radius Rj such that Jd(Rj) — Jh- First, this is 
not a useful operational definition, because in this case 
cos0 < — 0.5 Jd/ Jh — —0.5, is a static condition, which 
does not depend on the properties of the black hole or of 
the accretion disc. Second, matter contained within ra- 
dius Rj cannot transfer all its angular momentum to R w 
but only a fraction \/ ' R w / Rj Jd(Rj)- Therefore a more 
natural radial scale in the problem is the warp radius R w 
and in the following we will assume that Jd = Jd(Rw) 
(note that we share this choice with Lodato & Pringle 
2006). 

3. METHOD 

We explore the dependence of the alignment timescale 
in a Shakura-Sunyaev disc on: viscosity ^2/^1, black 
hole mass Mbh, misalignment angle, Eddington ratio, 
accreted mass m. Some articles on this subject use the 
solution of Collin-Souffrin & Dumont (1990), however in 
the view of the basic uncertainties of the problem we 
decided to use the less refined solution of Shakura & 
Sunyaev 1973. 

Assuming a Shakura-Sunyaev disc ("middle region"), 
the warp radius (in units of the Schwarzschild radius R$) 
can be expressed as: 
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where f Edd = Mc 2 /L Edd - 
We can then define the accretion timescale: 
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(where v\ = aH 2 flx was used), and the timescale for 
warp propagation: 
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The ratio of angular momenta of the disc at R Wl defining 
Md(Rw) = Mt acc (R w ), and of the black hole is: 
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4. RESULTS 

4.1. Single accretion episodes 

We first discuss the behavior of the disc+BH system 
during the alignment process. The scheme we adopt is 
as follows: 




Fig. 1. — Evolution of the misalignment angle between the angu- 
lar momentum vector of the outer accretion disc and the BH spin 
(top panel), and of the magnitude of the BH spin (bottom panel), 
due to accretion of aligned material (spin-up) or counter-aligned 
material (spin-down). The initial BH mass is 

Mbho = 5 x 1O 6 M , 

the initial spin a = 0.5, ui = v\, a = 0.1, and the accretion rate 
is /Edd = !■ The four curves show different initial misalignment 
angles (top to bottom: <j> = 3, 2, 1, 0.5 radians. ) 

1. for a BH with initial mass Mbho determine the 
initial conditions: warp radius, R w , timescale for 
warp propagation t w — R^/vi (y^ is chosen ei- 
ther coincident with v\, or V2 = v\/a 2 ), accre- 
tion timescale for material at the warp radius, 
iacc = R w /vx, angular momentum of the hole and 
of the disc at R w , Jh and Jd- 

2. using the King et al. (2005) condition for mis- 
alignment determine if the BH and the inner disc 
are aligned or counter-aligned (counter-aligned if 
cos0 < -0.5 Jd/Jh)- 

3. over timesteps At — t acc (R w ) compute the neces- 
sary quantities at the end of every step: increase in 
black hole mass due to accretion, new BH spin (fol- 
lowing Bardeen 1970; where the counter-alignment 
or alignment conditions are taken into considera- 
tion, i.e., BHs can be spun down or up), new Jh, 
new R w , new Jd, new angle between Jh and Jd 
(vectorial sum). In every timestep the disc within 
R w is consumed. 

Figures [Tj and [5] give examples of single accretion 
episodes for different initial angles between the angular 
momentum vector of the outer (not warped) portions of 
the accretion disc and the black hole spin. They show 
the evolution of the spin magnitude and inclination as 
computed for v 2 = v\ (Fig.l) and z^/^i = 1/a 2 (Fig. 2). 

As we can see, the alignment timescale is basically in- 
dependent of the misalignment angle. To modify signif- 
icantly the BH spin, one has to bring to R w an amount 
of angular momentum comparable to Jh- Therefore, if 
J h > Jd{R w ), then: 
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Since J h oc aM B K\/Rs, Jd{R w ) oc M d (R w )y/R^, and 
M d (R w ) = Mt acc (R w ), Equation E] gives (Rees 1978): 
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Denning the mass accreted during i a iign as m a ii gn = 
iaiignM, one gets: 
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Therefore a series of many randomly oriented accretion 
events with accreted mass m <C Tiaiign should result in 
black-hole's spin oscillating around zero. For the oppo- 
site case of m ^> realign the black hole will be spun-up 
to large positive spins; for m ~ Mbh the hole will be 
spun-up to d ~ 1. Let us notice finally that since it is 
reasonable to assume that m a ii g n "C Mbh the existence of 
AGN hosting black holes with d ~ — 1 is rather unlikely 
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Fig. 2. — As in Figured] but uq./v\ = 1/a 2 , a = 0.1. 



Our calculations show that it is difficult to avoid high 
spin for the most massive black holes. For large BH 
masses, the accretion timescale is very long, and con- 
sequently the warp radius and the angular momentum 
within the warp radius, J d , are large. If J d > 2 Jh, then 
the value |0.5 J d /Jh\ > 1 and the counter-alignment con- 
dition cannot be satisfied for any angle. This condition 
corresponds to: 
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If V2 = v\, Mbh, max is of order 10 7 — 10 s M Q for most 
sensible choices of a and /Edd- In this case the most 
massive black holes force accretion to occur from aligned 
discs, therefore causing a systematic spin-up. If the warp 
propagation is instead better described by v% = vi/a 2 , 
Mbh, max becomes exceedingly high and large accretion 



events can still act to spin down the black hole, provided 
m < m a i ign . 

The condition expressed in Equation [5] is true only if 
there is enough mass to fill the warp radius, that is if the 
total mass of the disc is larger than: 
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If the mass to be accreted by the BH in an episode is 
smaller than M(j jm i n , then J d <C Jh, and both alignment 
or counter-alignment can happen. 

4.2. Multi- accretion events 

We then run a series of simulations in which we ex- 
plore different parameters. We start with a small BH, 
Mbho = 1O 5 M , and have it grown by a series of ac- 
cretion episodes. The accreted mass m is randomly ex- 
tracted from only one of two different distributions: (1) 
a distribution flat in m, with m < O.IMbh, (2) a dis- 
tribution flat in m, with m < O.OIA/bh- The angle 
<p is extracted from a flat distribution < <f> < ir at 
the beginning of every accretion episode. Every sim- 
ulation is composed by a large number of accretion 
episodes, until one of the following conditions are met: 
M B h > 10 9 M Q or t tot > 10 10 years, that is the total 
simulation time (total time a BH accretes to grow from 
the initial Mbho = 10 5 M© mass to its final mass) is 
shorter than the age of the universe. During an episode 
where the BH accretes counter-aligned material, the BH 
is spun down. If the black hole is spun-down until its spin 
is zero, any subsequently accreted matter acts to spin the 
BH up again, although the direction of the spin axis is 
now reversed and aligned with the angular momentum 
of the disc. 
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Fig. 3. — Distribution of BH spins in different mass ranges. The 
accreted mass m is randomly extracted from a distribution flat in 
m, with m < O.IMbh- Initial mass Mbho = 10 5 Mq , initial spin 
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We run 100 simulations for every parameter sets choice, 
and we trace the spins at the end of every accretion 
episode, for all the accretion episodes in the simulations. 

We have explored a wide range of parameters, and we 
summarize here our findings. We have varied the accre- 
tion rate, from /Edd = 0.05 to /Edd = 1- If the accretion 
rate is low, the main caveat is that black holes do not 
reach high masses within the Hubble time, however, the 
efficiency of alignment is not strongly dependent on /Edd 
(see Equation [5]) . 

We have considered different black hole spins at birth, 
from a = 10 -3 to a — 0.9. After the BHs have changed 
their initial mass by about one order of magnitude, the 
distributions are indistinguishable from each other. Dur- 
ing the first few e-foldings, however, the spin distribution 
is peaked around the black hole spin at birth. 

We have also varied the viscosity parameter a (see sec- 
tion 2.1), and the relation between the viscosity charac- 
terizing the warp propagation {y^) with respect to the 
viscosity responsible for the transfer of the component of 
the angular momentum parallel to the spin of the disc 
(vi). When a is varied, but ^2/^1 is kept fixed, the dif- 
ferences between the spin distributions are not large. A 
smaller a skews the distribution towards higher spins (cf 
Equations [3 El). 
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Fig. 4. — As in Figure [3] but with m randomly extracted from 
a distribution flat in m, with m < O.OIMbh- 



One of the main parameters influencing the spin dis- 
tribution is the relation between 1/2 and v\. If vijv\ = 1, 
after the BHs have reached ttibh ~ 10 6 M Q , the spin dis- 
tribution is dominated by rapidly spinning black holes. 
Equation [5] also shows that the most massive black holes 
force accretion to occur from aligned discs, therefore 
causing a systematic spin-up, unless very small parcels of 
material are accreted at every single accretion episode. 
If the warp propagation is instead better described by 
a high V2 = Ui/ar, MBH.max becomes exceedingly high, 
and all sorts of accretion events can still act to spin down 
the black hole, provided m < m a ii gn . 

In fact, we confirm the results by MSL, that is that the 
main parameter governing the distribution of BH spins 



is the amount of material accreted in a single accretion 
episode. This result is clear from Figures H which refer 
to different choices for the distributions of m. Only if the 
mass accreted in one episode is smaller than m a ij gn , the 
distribution of black hole spins can remain flat. 

In the next section we discuss the likelihood of different 
m distributions in the light of evolutionary models for the 
BH population in a hierarchical cosmology. 

4.3. Merger driven accretion 

We first present an evolutionary track for BH spins, 
where a BH grows by a sequence of randomly oriented 
accretion episodes in a merger driven scenario. The BH 
mass evolutionary tracks are extracted from semianalyt- 
ical simulations of BH growth that have been shown to 
reproduce the evolution of the BH population as traced 
by the luminosity function of quasars (Marulli et al. 
2006, 2007; Volonteri, Salvaterra & Haardt 2006). We 
focus here on two specific tracks, for a putative BH in 
an "elliptical" (E) galaxy, and one in a putative "disc" 
(D) galaxy (Fig [5]). Here the morphological classifica- 
tion is purely based on the frequency of major mergers, 
i.e., mergers between comparable mass galaxy systems 
which are believed to contribute mainly to the spherical 
component of galaxies. A BH hosted in an "elliptical" 
galaxy should have experienced a major accretion event 
in connection with the last high-redshift major merger, 
which formed the elliptical galaxy as we see it now. Af- 
terwards, the galaxy (BH) has not grown in mass due 
to merger driven star formation (accretion). In the case 
of the BH hosted in a "disc" galaxy, a small number of 
minor mergers might have happened after the last major 
mergers. These minor mergers are believed to be respon- 
sible for re-building the galaxy disc. In conjunction with 
these minor mergers, a small infall of gas can produce a 
relatively minor accretion episode onto the BH as well. 




4 

z 



Fig. 5. — Growth of BHs in putative "elliptical" (upper panel) or 
"disc" (lower panel) galaxy ( only merger driven accretion events 
are considered). The initial spin a = 10 -3 , V2 = vi/ct 2 , a = 0.03, 
and the accretion rate is /Edd =0.1. 

Along the evolutionary tracks for our BHs, we trace the 
joint evolution of accretion onto the BH, the dynamics 
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vations suggest that single accretion events last ~ 10 5 
years in Seyfert galaxies, while the total activity lifetime 
(based on the fraction of disc galaxies that are Seyfert) is 
10 s - 10 9 years (e.g., Kharb et al. 2006; Ho et al. 1997). 
This suggests that accretion events are very small and 
very 'compact'. 

A type of random event which leads to short-lived ac- 
cretion episodes is the tidal disruption of stars. One 
expects discs formed by stellar debris to form with a 
random orientation. Stellar disruptions would therefore 
contribute to the spin-down of BHs. Let us consider the 
maximal influence that feeding via tidal disruption of 
stars can have on spinning down a BH. The number of 
tidal disruptions of solar type stars in an isothermal cusp 
per billion years can be written as: 



N* = 4 x 10 5 
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Fig. 6.— Spin of BHs at z 
"disc" (lower panel) galaxies. 



in "elliptical" (upper panel) or 



of the accretion disc, and the consequences on the spin. 
The scheme we adopt is similar to the one described in 
Section 4.1. During an episode where the BH accretes 
counter-aligned material, the BH is spun down, until the 
spin is zero, and subsequently any accreted matter acts 
to spin the BH up again, although the direction of the 
spin axis is now reversed and aligned with the angular 
momentum of the disc. 

In case "E" , the last accretion episode caused a large 
increase in the BH mass, following the major merger 
which created the elliptical itself (Hopkins & Hernquist 
2006). During this episode the spin increased signifi- 
cantly as well, up to very high values. Let us remind here 
that in the extreme event of a maximally-rotating hole 
spun down by retrograde accretion, the BH is braked af- 
ter its mass has increased by the factor ^3/2. Any mass 
accreted afterwards spins up the black hole, and if the 
final mass increase is by a factor 3, the BH will end up 
maximally-rotating again. 

Let us now consider a "D"-type evolution. A BH would 
experience a series of small accretion episodes (triggered 
possibly by minor mergers), extending for a longer period 
of times. If these episodes are uncorrclated, that is if 
the inflow during a given episode is not aligned with the 
orientation of the spin of the BH, the randomization of 
the angle </> over the (few) accretion episodes tends to 
spin down the BH. 

We run a statistical sample of "E" and "D" track, 
for BHs hosted in large (i.e., Andromeda-size systems) 
galaxies. We find that, if only the merger driven evo- 
lution is taken into account, BHs in ''elliptical" galaxies 
are left with large spins. BHs in "disc" galaxies have, on 
average, slightly lower spins, however the distribution is 
still peaked at large values (Figure [5]). 

4.4. Short-lived accretion events 

Even minor mergers tend to trigger inflows of matter 
which arc too large to lead to the series of short lived 
accretion events necessary to leave BHs with small spins 
(cfr. the discussion in MSL). Moreover, several obser- 



Assuming that BH masses scale with the velocity dis- 
persion, er, of the galaxy (we adopt here the Tremaine et 
al. 2002 scaling), we can derive the relative mass increase 
for a BH in 1 billion years: 



A/bh 



0.37 



/ A/] 



BH 



V 10 6 Mr. 



-9/8 



(12) 



The maximal level of spin down would occur assuming 
that all the tidal disruption events form countcrrotat- 
ing discs, leading to retrograde accretion (note that the 
mass of the debris disc is much smaller than M<j im i n , cfr. 
Eq. [TQl so that counter-alignment is allowed for any BH 
mass). Eq. [T2l shows that a small (say 10 6 M Q ) BH start- 
ing at a — 0.998 would be spun down completely, on the 
other hand the spin of a larger (say 10 7 M Q ) BH would 
not be changed drastically. 

Early type discs typically host faint bulges character- 
ized by steep density cusps, both inside (Bahcall & Wolf 
1976; Merritt & Szell 2006) and outside (Faber et al. 
1997) the sphere of influence of the BH. In this envi- 
ronment, the rate of stars which are tidally disrupted by 
BHs (Hills 1975; Rees 1988) less massive than 10 8 M Q 6 is 
non negligible (Milosavljevic et al. 2006). Since in ellip- 
tical galaxies the central relaxation timescale is typically 
longer than the Hubble time, and the central density pro- 
file often displays a shallow core, tidal disruption of stars 
is unlikely to play a dominant role. 

An additional feeding mechanism might be at work in 
gas-rich galaxies with active star formation. Compact 
self-gravitating cores of molecular clouds (MC) can oc- 
casionally reach subparsec regions, and may do it with 
random directions provided that the galactic disc is much 
thicker than the spatial scale of the BH gravity dom- 
ination region (Shlosman, private communication). Al- 
though the rate of such events is uncertain, we can adopt 
the estimates of Kharb et al. (2006), and assume that 
about 10 4 of such events happen. We can further assume 
a lognormal distribution for the mass function of MC 
close to galaxy centers (based on the Milky Way case, 
e.g., Perets, Hopman & Alexander 2006). We do not 
distinguish here giant MC and clumps, and, for illustra- 

6 For black hole masses ^2x 10 8 Mg the Schwarzschild radius 
exceeds the tidal disruption radius for main-sequence stars. 
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2xlO a 4x10 s 6xl0 8 8x10 s 
total acc. time (yr) 

Fig. 7. — Evolution of BH spins due to accretion of molecular 
clouds cores. We assume a lognormal distribution for the mass 
function of molecular clouds (peaked at Io%{Mmc / Mq) = 4, with 
a dispersion of 0.75. The initial spin of the BHs is 0.998. Upper 
curve: the initial BH mass is 1O 7 M0, lower curve: the initial BH 
mass is 10 s M Q 

tive purpose we assume a single lognormal distribution 
peaked at log(MMc/ M Q ) = 4, with a dispersion of 0.75. 

Fig [7] shows the possible effect that accretion of molec- 
ular clouds can have on spinning BHs. The result is, on 
the whole, similar to that produced by minor mergers 
of black holes (Hughes & Blandford 2003), that is a spin 
down in a random walk fashion. The larger the BH mass, 
the more effective the spin down. 

In a gas-poor elliptical galaxy, however, substantial 
populations of molecular clouds are lacking (e.g. Sage 
et al. 2007), thus hampering the latter mechanism for 
short lived accretion events proposed. 

5. DISCUSSION AND CONCLUSIONS 

We have investigated the evolution of BH spins driven 
by accretion from discs with angular momentum vectors 
that can be misaligned with respect to the spin axis. 
We have assumed that accretion discs can be described 
by Shakura-Sunyaev a-discs, and that when the angu- 
lar momentum of the accretion disc is not aligned with 
the spin of the BHs, the disc itself is warped. The 
inner portions of the discs experience Lense-Thirring 
torque, which tends to align the inner parts of the disc. 
The timescale of the Lense-Thirring precession increases 
faster with distance from a BH than the timescale of 
warp propagation, and they equate at the so-called "warp 
radius" , where a transition occurs from alignment to 
misalignment. King et al. (2005) pointed out that for 
highly misaligned discs, counter-alignment, rather than 
alignment can occur. The co- or counter-alignment of 
the accretion discs has important consequences on the 
spin of BHs. A black hole accreting from prograde or- 
bits (i.e., alignment case) is spun up by the coupling 
between the angular momentum of the infalling mate- 
rial and its spin (Bardeen 1970). If, instead, an ini- 
tially spinning hole accretes from retrograde orbits (i.e., 
counter- alignment case), it is spun down. An initially 
non-rotating BH gets spun up to a maximally-rotating 



state (d = 1) after reaching the mass Mbh = VoMbho- 
A maximally- rotating hole (d = 1) gets spun down by 
retrograde accretion to d = after reaching the mass 
M BH = T/3/2M Bm . A 180° flip of the spin of an 
extreme-Kerr hole will occur after Mbh = 3Mbho- 

It is therefore necessary that accretion episodes in- 
crease the mass of a BH by less than Mbh = VGMbho, 
in order to keep the spin at low values, if accretion 
preferentially occurs from prograde orbits. Natarajan 
& Pringle (1998) suggested that accretion indeed oc- 
curs from aligned discs (i.e. prograde orbits), as the 
timescale for disc alignment is much shorter than the 
timescale of the BH mass growth by Mbh = V^Mbho- 
King et al. (2005) suggested, however, that when the 
initial misalignment angle is large and m is sufficiently 
small, counter-alignment, rather than alignment, occurs 
and BHs can be spun down in a large fraction of the 
accretion episodes. 

We have quantified here the likelihood of counter- 
alignment and spin down as claimed by King et al. 
(2005). We identify two main parameters influencing 
the distribution of BH spins: the distribution of the ac- 
creted mass, m, with respect to the mass of the BH, and 
the relation between 1/2 and v\, where V2 is the viscos- 
ity characterizing the warp propagation, and v\, which 
is responsible for the transfer of the component of the 
angular momentum parallel to the spin of the disc. V2 
can in principle differ from v\. 

If the accreted mass, m, is much smaller than the mass 
of the BH (e.g., m < 0.01 Mbh), the distribution of black 
hole spins is flat, as the timescale for spin overflipping 
due to the Bardeen-Petterson effect is longer than the 
timescale to accrete the whole m. If instead m ~ Mbh 5 
BHs can align with the angular momentum of the accre- 
tion disc, and accrete enough mass to be spun up. In 
this case the distribution of BH spins is dominated by 
rapidly rotating systems. 

Understanding if the description of the warp propaga- 
tion is correctly described by a different viscosity with 
respect to the one responsible for the radial propagation 
of the angular momentum is beyond the scope of this 
work. We have therefore explored a wide range of pos- 
sible viscosities, and we simply report here our results. 
If V2jv\ = 1 the timescale for alignment is short, and 
the spin of a BH increases rapidly. If the warp propaga- 
tion is instead better described by a high vi = v\/a 2 , a 
substantial fraction of black holes of all masses can have 
small spins, provided m <C Mbho- 

However, both semi-analytical models of the cosmic 
BH evolution (Volonteri et al. 2005) and simulations of 
merger driven accretion (di Matteo et al. 2005) show that 
most BHs increase their mass by an amount 3> m a n gni 
if the evolution of the LF of quasars is kept as a con- 
straint. These high m values are likely characteristic of 
the most luminous quasars and most massive black holes 
- especially at high redshift. We expect therefore that 
bright quasars at z > 3 have large spins, in contrast with 
the suggestion of King & Pringle (2006). High spins in 
bright quasars are also indicated by the high radiative 
efficiency of quasars, as deduced from observations ap- 
plying the Soltan argument (Soltan 1982; Wang et al. 
2006 and references therein). 

If the mass of a BH need to reach 10 9 M by z — 3, or 



even more strikingly, by z — 6, so that they can represent 
the engines of quasars with luminosity L > 10 46 ergs _1 , 
BHs need to grow from typical seed masses (e.g. Madau 
& Rees 2001, Koushiappas et al. 2003, Begelman, Volon- 
teri & Rees 2006, Lodato & Natarajan 2006) by at least 
3-4 orders of magnitude in 10 s — 10 9 years. The necessity 
of long and continuous accretion episodes implies there- 
fore that, for these BHs, m ^> m a i ign . 

Smaller BHs, powering low luminosity Active Galactic 
Nuclei, can instead grow by accreting smaller packets of 
material, such as tidally disrupted stars (for BHs with 
mass < 2 x 10 6 M Q , Milosavljevic et al. 2006), or pos- 
sibly molecular clouds (Hopkins & Hernquist 2006). For 
these black holes the spin distribution is more probably 
flat, or skewed towards low values. This latter result is 
in agreement with Sikora, Stawarz & Lasota (2006), who 
find that disc galaxies tend to be weaker radio sources 
with respect to elliptical hosts. In the hierarchical frame- 
work we might expect that the BH hosted by an ellipti- 
cal galaxy had, as last major accretion episode, a large 
increase in its mass following the major merger which 
created the elliptical itself (Hopkins & Hernquist 2006). 
During this episode the spin increased significantly as 
well, possibly up to very high values. Subsequently the 
black hole might have grown by swallowing the occasional 
molecular cloud, or by tidally disrupting stars. If the to- 
tal contribution of these random episodes represents a 
small fraction of the BH mass, the spin is, however, kept 
at high values. 

Black holes in spiral galaxies, on the other hand, prob- 
ably had their last major merger (i.e., last major accre- 
tion episode), if any, at high redshift, so that enough 
time elapsed for the galaxy disc to reform. Most of the 
latest growth of the BH should have happened through 



minor events, which have likely contributed to the BH 
spin down. 

Our results are supported also by the recent finding 
by Capetti & Balmaverde (2006; 2007) that radio bi- 
modality correlates with bimodality of stellar brightness 
profiles in galactic nuclei. The inner regions of radio loud 
galaxies display shallow cores (star deficient). Cores, 
in turn, are preferentially reside in giant ellipticals (see 
Lauer et al. 2007 and references therein). Radio quiet 
galaxies, including nearby low luminosity Seyferts, have 
instead power-law (cuspy) brightness profiles and prefer- 
entially reside in SO and spiral galaxies. 

Hence, noting that core nuclei result from merging 
BHs following galaxy mergers (Ebisuzaki et al. 1991, 
Milosavljevic & Merritt 2001, Milosavljevic et al. 2002, 
Ravindranath et al. 2002, Volonteri et al. 2003), Bal- 
maverde & Capetti's discovery is consistent with our 
conjecture that spin bimodality is determined by diverse 
evolutionary tracks of BH spins in disc galaxies (ran- 
dom small mass accretion events) and giant elliptical 
galaxies (massive accretion events which follow galaxy 
mergers). Both tidal disruption of stars, and accretion 
of gaseous clouds is unlikely in shallow, stellar domi- 
nated galaxy cores. Therefore it is conceivable that the 
observed morphology-related bimodality of AGN radio- 
loudncss results from bimodality of central black-holes 
spin distribution. 
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